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ABSTRACT 

We present a catalog of 1750 massive stars in the Large Magellanic Cloud, 
with accurate spectral types compiled from the literature, and a photometric cat- 
alog for a subset of 1268 of these stars, with the goal of exploring their infrared 
properties. The photometric catalog consists of stars with infrared counterparts 
in the Spitzer SAGE survey database, for which we present uniform photometry 
from 0.3 - 24 //m in the UBVIJHK S +IRAC+MIPS24: bands. The resulting 
infrared color-magnitude diagrams illustrate that the supergiant B[e], red super- 
giant and luminous blue variable (LBV) stars are among the brightest infrared 
point sources in the Large Magellanic Cloud, due to their intrinsic brightness, 
and at longer wavelengths, due to dust. We detect infrared excesses due to free- 
free emission among ~ 900 OB stars, which correlate with luminosity class. We 
confirm the presence of dust around 10 supergiant B[e] stars, finding the shape 
of their spectral energy distributions (SEDs) to be very similar, in contrast to 
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the variety of SED shapes among the spectrally variable LBVs. The similar lu- 
minosities of B[e] supergiants (logL/L > 4) and the rare, dusty progenitors 
of the new class of optical transients (e.g. SN 2008S and NGC 300 OT), plus 
the fact that dust is present in both types of objects, suggests a common origin 
for them. We find the infrared colors for Wolf-Rayet stars to be independent of 
spectral type and their SEDs to be flatter than what models predict. The results 
of this study provide the first comprehensive roadmap for interpreting luminous, 
massive, resolved stellar populations in nearby galaxies at infrared wavelengths. 

Subject headings: infrared: stars- stars: early-type- stars: Wolf-Rayet- stars: 
emission-line, Be- galaxies: individual (LMC)- catalogs 



Introduction 



The Spitzer Sp ace Telescope Legacy Survey called "Surveying the Agents of a Galaxy's 
Evolution" (SAGE; iMeixner et al.l 120061 ) provides an unprecedented opportunity to inves- 
tigate the infrared properties of a large number of massive stars (>8MjJ in the Large 
Magellanic Cloud ( LMC). The SAGE infrared point source catalog has so far enabled stud- 



ies of ev olved stars (Blum et al. 2006 : Hora et aL 2008 : Srinivasan et al.l2009l ). young stellar 



objects (jWhitney et al.l 120081 ) . and variables (jVijh et al.l 120091 ); however, the hot massive 
star population remains unexplored. The radiation from hot massive stars at the mid- 



infrared wavelengths probe d by the IRAC (3.6, 4.5, 5.8, 8.0 /xm; iFazio et al.l 12004 ) and 
MIPS (24, 70, and 160 /im; iRieke et al.ll2004l ) cameras onboard Spitzer consists of: thermal 
blackbody emission modified by the atmospheric opacity, bound-free and free-free emission 
("Bremsstrahlung"), which depend on the density and velocity structure of the stellar wind, 
and excess emission if circumstellar dust (or a cool companion) exists. 



Panagia fc Fellil (119751 ) and IWright fc Barlow! (119751 ) were the first to calculate the in 



frared and radio free-free emission from ionized envelopes around hot massive stars, as a 
function of the mass- loss rate (M), terminal velocity (t>oo) of the wind and temperature 
(T) for the cases of optically thin and thick winds to explain the radio emission from 
hot stars undergoing mass-loss. They found that the flux as a function of frequency (v) 
and distance (d) for an optically thick, spherically symmetric wind at infrared and radio 
wavelengths scales as: F ut hi C k oc M f t ^, 4//3 ^°' 6 d _2 T ' 1 , while for an optically thin wind: 



vlhin 



oc M Vo^R- 1 T 



-1.35, 



21 B u , where R is the photospheric radius, v Q the initial wind 



These stars are massive enough to fuse iron in their cores; they end their lives as core-collapse supernovae. 
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velocity at the surfac e of the star, and B y is the Planck function. These calculations, follow- 
ing the shell model of iGehrz et al.l (119741 ). motivated further infrared and radio observations 
of Galactic O B stars. The initia l studies produced controversial measureme nts of infrared 



excesses see 



Sneden et al. 1978; Castor fc Simon 1983 



Abbott et al.l 19841). whereas un- 



ambig uous excesses were detected in the radio (e.g. lAbbott et al.lll98ll ). iLeitherer fc Woli 
( 119841 ) detected infrared excesses in a large number of early-type stars and attributed the pre- 
vious non-detections to their different treatment of intrinsic stellar colors and interstellar red- 



dening. iBertout et al.l (119851 ) determined the velocity law for this large sample, pointing out 
that excesses can only be measured reliably in cases of suitably slow (low Voo) and/or dense 
(high M) winds. In principle, mass-loss rates can be determined from radio observations , 
which probe the optic ally thick wind, e.g . as done for OB-ty pe stars (lAbbott et ajj]l 98l[) 



and Wolf-Rayet stars (Barlow et al.lll98ll : lAbbott et al.lll986f ). However, Moffat fc Robert 



( 119941 ) presented observational evidence for clumping in stellar winds, which yields mass-loss 
rates that are too high when unaccounted for. While clumping, when not accounted for, can 
undoubtedly lead to an overestimation of the mass loss, the radio observations seem to be 
least affected by this problem (see e.g., iPuls et al.l 120061 ) . possibly because the wind is less 
clumped at larger distances from the star. On the other hand, it has been shown that in 
many c ases a substantial fraction of the ra dio flux from early-type stars is of non-thermal 
origin (lAbbott et al.lll984l ; |Pe Becker! 120071 ). thus complicating the interpretation of radio 
observations. 

Until now, work on infrared excesses in hot stars has been limited by (i) the accuracy 
of ground based observations at wavelengths greater than 5 /im, where infrared excess due 
to free-free or dust emission becomes significant, and (ii) the systematic errors associated 
with disentangling reddening and distances for Galactic stars, which can be substantial. The 
advent of Spitzer, therefore, provides an opportunity to readdress and quantify the infrared 
excesses in hot massive stars at half-solar metallicity. For the first time, a large number of 
stars in the LMC have been observed uniformly out to 24 /im. The low foreground reddening 
of this sample of stars, located at the same distance, circumvents the problems encountered 
with extincted Galactic early-type stars, which are observed against a bright infrared back- 
ground and distributed at various distances throughout the Galactic plane. With the goal of 
characterizing their infrared properties, we have compiled a catalog of luminous massive stars 
with known spectral types in the LMC and extracted infrared photometry for these stars 
from the SAGE database. Our results not only serve as a valuable starting point for more 
detailed work on specific types of massive stars, but also provide a roadmap for interpreting 
existing and future infrared observations of the most luminous stars in nearby galaxies. 

Besides the LM C, Spitz er has imaged the followin g Local Group ga l axies: the SMC 



dBolatto et al.l 120071 ). M31 dMould et al.l 120081 ). M33 (IVerlev et all 120071 : McQuinn et al. 
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NGC 6822 flCannon et al. 



20071: iThompson et al.ll2008h. W LM fl Jackson et alJl2007ah . IC 1613 fl Jackson et al.ll2007bh 



2006 ), and the dw arf irregulars: Phoenix, LGS 3, DDO 210, Leo A, 



Pegasus, and Sextans A (IBoyer et al.l 120091 ) . Only a small fraction of the massive stars in 



these galaxies have spectroscopic classifications, rendering identifications of objects based on 
their infrared colors uncertain. Our results, which are based on stars with accurate spectral 
classifications, will help interpret these photometric observations. An advantage of study- 
ing stars in the LMC, rather than in more distant Local Group galaxies, is that blending 
problems are minimized. For example, the angular resolution at the IRAC 3 .6 £tm band 
is l."7, which corresponds to 0.4 pc in the LMC (at a distance of 48 .1 kpc; iMacri et al. 



20061 ). and 8 pc in M33 (at a distance of 964 kpc; iBonanos et al.l 120061 ) . while typical OB 
associations have sizes of tens of parsecs (but half lig ht radii of several parsecs), and young 
massive clusters can have sizes smaller than a parsec (IClark et al.l 120051 ) . 



The paper is organized as follows: §2] describes our spectroscopic and photometric cata- 
logs of massive stars in the LMC, £J3] presents the resulting color-magnitude diagrams, §Hthe 
two-color diagrams, and $5] the infrared excesses detected in OB stars. Sections 9SI 9ZI 
§9] describe the infrared properties of Wolf-Rayet stars, luminous blue variables, supergiant 
B[e] stars, and red supergiants, respectively, and CTU1 summarizes our results. 



2. Catalog of Massive Stars in the LMC 

We first compiled a catalog of massive stars with known spectral types in the LMC from 
the literature. We then cross-matched the stars in the SAGE database, after incorporating 
optical and near-infrared photometry from recent surveys of the LMC. The resulting photo- 
metric catalog was used to study the infrared properties of the stars. In §2. II we describe the 
spectral type catalog compiled from the literature, in §2.21 the existing optical, near-infrared 
and mid-infrared surveys of the LMC that were included in the SAGE database, and in §2.31 
the cross-matching procedure and the resulting photometric catalog. 



2.1. Spectral Type Catalog 



The larges t existing catalogs of O B stars in the LMC were compiled by lSanduleak! (jl970j, 
1273 stars) and lRousseau et al.l (119781 . 1822 stars); however the accuracies of the coordinates 
(arcminute precision) and spectral types (from objective prism spectra) are not sufficient for a 
comparison with the SAGE database. Despite the monumental effort of Brian Skiff to update 
coordinates of these stars by manually identifying each one from the original fmdercharts (see 
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Sanduleakll2008l ) . improved spectral types require new spectroscopy and acquiring them is an 
even more difficult task. We therefore compiled a new catalog of massive stars from the lit- 
erature, by targeting OB stars having both accurate coordinates and spectral classifications, 
making some of these available for the first time in electronic format. The literature search 
resulted in 1750 entries. We do not claim completeness; however, we have targeted studies 
of hot massive stars (defined by their spectral types) in individual clusters and OB associa- 
tions, as wel l as st udies of particular types of massive stars . The larges t studies incl uded are: 
Conti et al. (1986, new spectral ty pes for 191 OB stars), iFitzpatrickl (jl988l, 19~9~ll. study of 
supergiants , Massey et al. ( 1995 . new spectral types for 179 OB stars ), Massev et al. ( 2000 



study of 19 OB associations), the VLT-FLAMES survey flEvans et alJl20oJ 2 38 stars in Nil 



and N GC 2004), the studies o f the 3 Doradus region by ISchild fc Tested <ll992n. IParker 
( 1993 ) and IWalborn fc Blades! (119971 ). and the Wolf-Rayet catalog from iBreysacher et al. 
(119991 ). Note, we have omitted st ars from the crowded R136 clu ster. We have added the 



6 known lumin ous blue variables ([Humphreys fc Davidson 



2005 



19941). 10 supergian t B[e] stars 



Raguzova fc Popov! 120051 ) and early 



(jZickgrafj 120061 ). 20 Be/X-ray binaries (ILiu et al. 

f 

(1l979l ). The completeness of the catalog depends on the spectral type, e.g. it is ~ 3% for 



type eclipsi ng and spectroscopic bi naries. For comparison, we also included 14 7 red super- 
giants from iMassey fc Olsenl (120031 ) , some of which were originally classified by Humphreys 



the unevolved O stars in our catalog (out o f an estimated total of 6100 une volved stars with 



mass es > 20 M Q in the LMC; iMasseyl 120091 ). while the Wolf-Rayet catalog (IBreysacher et al. 
19991 ) is thought to be close to complete. 



Table [T] presents our catalog of 1750 massive stars, listing: the star name, coordinates 
in degrees (J2000.0), the reference and corresponding spectral classification, along with com- 
ments. The names of the stars are taken from the corresponding reference and are sorted 
by right ascension. The spectral classifications in the catalog are typically accurate to one 
spectral type and one luminosity class. We retained about a dozen stars from the above 
studies with only approximate spectral types (e.g. "early B"). For double entries, we in- 
cluded the best spectral classification available, usually corresponding to the most recent 
reference. We updated coordinates for all stars with names from the Sanduleak, Brunet, and 
Parker catalogs from Brian Skiff's lists^l Our catalog contains 427 stars from the Sanduleak 
catalog, 81 from the Brunet catalog and 148 from the Parker catalog. 



2 We have corr ected the typographical error in the name of the B7 Ia- 
Fitzpatrickl (jl99lh to the correct name: Sk -69° 143. 



star Sk -67° 143 given by 



3 ftp://ftp. lowell.edu/pub/bas/starcats/ 
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2.2. Optical and Infrared Surveys of the LMC 



Several large optical and infrared photometric catalogs of the LMC have recently become 
available, enabling us to obtain accurate photometry for its massive star population in the 
wavelength region 0.3 — 8 iim and in some cases up to 24 fim. The optical surveys are: 
the UBVR catalog of iMasseyl (120021 ) with 180,000 bright stars in an area covering 14.5 



squar e degrees, the UBVI Magellanic Clouds Photometric Survey (MCPS; IZaritsky et al. 



20041 ) including 24 million stars in the central 64 square degress, and the OG LE HI catalog 
conta ining VI photometry of 35 million stars covering 40 square degrees Judalski et al. 



20081 ). The angular resolution for each survey is ~ 2. "6 for the catalog of IMasseyl (120021 ). 
~ l."5 for MCPS, and ~ l."2 for OGLE III. 

The exist ing near-infrared pho tometric catalogs include: the Two Micron All Sky Sur- 
vey (2MASS: Iskrutskie et alJ bood. extended by 6X2MASS) and the targeted IRSF survey 
(IKato et al.l 120071 ). which contains 14.8 million sources in the central 40 square degrees of 
the LMC. 2MASS has a pixel scale of 2."0 x 2."0, an average seeing of 2. "5 and limiting 
magnitudes of J = 15.8, H = 15.1 and K s = 14.3. IRSF has a pixel scale of 0."45 pixel -1 , 
average seeing of l."3, l."2, in the JHK S bands, respectively, and limiting magnitudes of 
J = 18.8, H = 17.8 and K s = 16.6. In the mid-infrared, the Spitzer SAGE survey uniformly 
imaged the central 7° x 7° of the LMC in th e IRAC and MIPS bands on two epochs in 2005, 
separated by 3 months (IMeixner et al.ll2006l ). The survey has recently produced a combined 
mosaic catalog of 6.4 million sources. IRAC, with a pixel scale of l."2 pixel -1 , yields an 
angular resolution of l."7 — 2."0 and MIPS at 24 jum has a resolution of 6". 

Given the variation in the depth, resolution and spatial coverage of these surveys, we 
included the available photometry for the massive stars in our catalog from all the MCPS 



OGLE III, 2 MASS, IRSF and SAGE catalogs. MCPS has incorporated the catalog of lMassev 



( 120021 ) for bright stars common to both catalogs. We note that due to problems with zero- 
point, calibration and PSF fitting of bright stars, MCPS photometry of bright stars should 
be used with caution. The authors warn that s tars "brighter than 13.5 mag in B or V are 
prone to substantial photometric uncertainty" (IZaritsky et al.l 120041 ) . Photometry of higher 
accuracy, particularly in the optical, exists in the literature for many of the stars in our 
catalog; however, it was not included in favor of uniformity. 
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2.3. Photometric Catalog 

2.3.1. Matching Procedure 

We used a near-final pre-release version of the SAGE catalog to obtain near- to mid- 
infrared photometry for our catalog of 1750 massive stars. The SAGE catalog includes the 
IRAC mosaic photometry source list, extracted from the combined epoch 1 and epoch 2 
IRAC images and merged together w ith 2MASS and the 2MASS 6X Deep Point Source 



catalog (6X2MASS; ICutri et alj|200j ), and the MIPS 24 /zm epoch 1 catalog. The IRAC 
and MIPS 24 /jm catalogs were constructed from the IRAC and MIPS 24 /im source lists 
after applying stringent selection criteria and are highly reliable. These catalogs are subsets 
of more complete, but le ss reliable arch i ves. More details on the data processing and data 
products can be found in Meixner et al. ( 2006 ) and the SAGE Data Description Document^. 



We extracted infrared counterparts to the massive stars in our list from the SAGE 
IRAC mosaic photometry catalog by performing a conservative neighbor search with a 1" 
search radius and selecting the closest match for each source. This procedure yielded mid- 
infrared counterparts for 1316 of the 1750 sources. The IRAC mosaic photometry catalog, 
MIPS 24 /im catalog, IRSF, MCPS, and OGLE III catalogs were crossed-matched in the 
SAGE database to provide photometry for sources over a wavelength range from 0.3 — 
24 /im. We used this "universal catalog" to extract multi-wavelength photometry for IRAC 
sources matched to the massive stars. Specifically, for IRAC sources with one or more 
matches in other catalogs (all but 5), we only considered the closest matches between sources 
from any two available catalogs (IRAC-MIPS24, IRAC-MCPS, MIPS24-MCPS, IRAC- 
IRSF, MIPS24-IRSF, IRAC-OGLE III, MIPS24-OGLE III), with distances between the 
matched sources of <1". For example, for a match between the IRAC, MIPS 24 /mi and 
MCPS catalogs (IRAC-MIPS24-MCPS), we applied these constrains on the IRAC-MIPS24, 
MIPS24-MCPS, and IRAC-MCPS matches. These stringent criteria were used to ensure 
that sources from individual catalogs for each multi-catalog match refer to the same star; 
they reduced the source list to 1262 sources. Six additional massive stars had matches within 
1" in the MIPS 24 /im Epoch 1 Catalog, but not in the IRAC Mosaic Photometry Catalog 
(IRACC). We supplemented the missing IRAC data: two sources had counterparts in the 
IRAC Mosaic Photometry Archive (IRACA) and 4 MIPS 24 /im sources had counterparts in 
the IRAC Epoch 1 Catalog (IRACCEP1). These IRAC sources were the closest matches to 
the MIPS 24 /im sources within 1". Five out of 6 sources also have matches within 1" in the 
MCPS catalog. Table [2] shows the breakdown of the matched stars to the catalogs, such that 



4 http://irsa.ipac.caltcch.cdu/data/SPITZER/SAGE/doc/ 
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5 stars were only matched to the IRACC, 88 only to the IRACC+IRSF catalogs etc. The 
above requirements reduced the source list to 1268 sources^. We supplemented the missing 
IRAC data for 6 sources: two sources had counterparts in the IRAC mosaic photometry 
archive and 4 MIPS 24 /im sources had counterparts in the IRAC epoch 1 catalog. We defer 
the discussion of misidentifications and blending to §9j 



2.3.2. Catalog Description 

Table [3] presents our final matched catalog of 1268 stars, with the star name, IRAC 
designation, UBVIJHK S +IR,AC+MIPS24 photometry and errors, reference paper, corre- 
sponding spectral classification and comments, sorted by increasing right ascension. Overall, 
the photometry is presented in order of shortest to longest wavelength. The 17 columns of 
photometry are presented in the following order: UBVI from MCPS, VI form OGLE III, 
JHK S from 2MASS, JHK S from IRSF, IRAC 3.6, 4.5, 5.8, 8.0 /im and MIPS 24 /an. A 
column with the associated error follows each measurement, except for the OGLE III VI 
photometry. Henceforth, JHK S magnitudes refer to 2MASS photometry, whereas IRSF 
photometry is denoted by a subs cript, e.g. Jirsf- All magnitudes are calibrated relative to 



Vega (e.g. see iReach et al.l 120051 . for IRAC). In Table HI we summarize the characteristics 
of each filter: effective wavelength A e ff, zero magnitude flux (in Jy), angular resolution, and 
the number of detected stars in each filter. 

The spatial distribution of our 1268 matched sources is shown in Figure [Q, overlayed 
onto the 8 /im image of the LMC. We find that the massive star population traces the spiral 
features of the 8 /im emission, which maps the surface density of the interstellar medium. 
Despite the fact that most studies we included from the literature targeted individual clusters, 
the spatial distribution of OB stars in our catalog is fairly uniform, with the exception of the 
Nil, NGC 2004 and 30 Doradus regions, which have been subjects of several spectroscopic 
studies. 



3. Color— Magnitude Diagrams 



We divide the matched stars into 9 categories according to their spectral types: O stars, 
early (B0-B2.5) and late (B3-B9) B stars (the latter have supergiant or giant luminosity 



Preliminary estimates indicate that matches against the final catalog will not increase the number of 
matched stars by more than 1%. 
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classifications, except for 2 B3V stars), spectral type A, F and G (AFG) supergiants, K and 
M red supergiants (RSG), Wolf-Rayet stars (WR), supergiant B[e] stars (sgB[e]), confirmed 
luminous blue variables (LBV) and Be/X-ray binaries. In Figures |5] and [31 we present 
infrared [3.6] vs. [3.6] — [4.5] and J — [3.6] color-magnitude diagrams (CMDs) for all the 
stars in the catalog, identifying stars in these 9 categories. The conversion to absolute 



magn itudes in all CMDs is based on a true LMC distance modulus of 18.41 mag (jMacri et al. 



20061 ). The locations of all the SAGE catalog detections are represented by the grey two 
dimensional histogram (Hess diagram). The red giants form the clump at [3.6] > 15 mag, 
while the vertical blue extension contains late-type LMC and foreground stars (free-free 
emission causes the OB stars to have redder colors). The asymptotic giant branch (AGB) 
stars are located at [3.6] ~ 10 mag, J — [3.6] ~ 2.5 mag. Immediately, one notices that the 
sgB[e], RSG, and LBVs stars are among the brightest stars at 3.6 /iin and occupy distinct 
regions in the diagrams. Most of the O and B stars are located along a vertical line at 
[3.6] — [4.5] ~ 0, as expected. The RSG have "blue" co lors because of the depression at 



[4.5] due to the CO band (see e.g. IVerhoelst et al.ll2009i ). The WR stars are on average 



0.3 mag redder than the OB stars because of their higher wind densities and, therefore, 
stronger free-free emission, while the sgB[e] stars are 0.6-0.7 mag redder and have similar 
brightnesses to the RSG. The Be/X-ray binaries are located among the "red" early-B stars. 
The LBVs have similar colors to the WRs, but are brighter. The late-B stars are brighter 
than the early-B stars because most of the former are luminous supergiants. The brightest 
of the AFG supergiants is the (only) G-type supergiant Sk — 69°30. A reddening vector for 
E(B — V) — 0.2 mag is shown in Figure [3] to illustrate the small effect of reddening, which 
decreases at longer wavelengths. 

In Figures H] and [51 we present [8.0] and [24] vs. [8.0] — [24] CMDs, respectively. At these 
wavelengths, only stars with dust are detected, as the sensitivity of Spitzer drops sharply, 
while the flux from hot stellar photospheres also decreases. Therefore only the dusty WR, 
LBVs, sgB[e], and RSG are detected, as well as certain OB stars and late B and AFG 
supergiants, which are likely associated with hot spots in the interstellar medium, and may 
well be interacting with it. See ^JSl El El and M for a more detailed description of their 
colors. Stars with cool dust are brightest at 24 /im, while emission from warmer dust peaks 
at 8 /mi. In these CMDs, the locations of the massive stars overlap with those of AGB 
stars, young stellar objects and planetary nebulae. Infrared colors alone are not sufficient to 
distinguish among these very different types of stars, illustrating the important diagnostic 
value of wide-baseline photometry. 



-10- 



4. Two— Color Diagrams 

In Figures El [7] and [8] we present two-color diagrams (TCDs) using the near and mid- 
infrared photometry from our matched catalog. We label stars according to their spectral 
types and overplot all the SAGE detections in grey as Hess diagrams. We find that the 
majority of the OB stars have colors near zero in all bands, as expected, while there is an 
almost continuous extension of OB star "outliers" toward the WR stars, due to thick winds 
or the presence of disks. Among these outliers are the Be/X-ray binaries and other emission 
line stars. The most conspicuous group of stars in all TCDs are the B[e] supergiants, which 
have large excesses, of ~ 4 mag in the K — [8.0] color. The RSG occupy a distinct position 
in the CMDs, because of their cooler temperatures, but are found to have a range of mid- 
infrared colors due to the amount of dust they contain (see §9]). We have overplotted 6 
simple theoretical models to guide the interpretation of the stars in these diagrams. These 
models are: 

(i) Blackbodies (BBs) at 30,000, 10,000, 5,000, 3,000, 1,700, 1,100, 800 and 500 K. BBs 
are a good approximation of stellar emission in the infrared for most stellar atmosphere 
spectral energy distributions (SEDs), with the possible exception of red supergiants, which 
have spectra dominated by emission lines. 

(ii) A power law model F u cx z/*, for a ranging from —1.5 to 2 in steps of 0.5. In the 
infrared, most stellar SEDs are c lose to the Rayle i gh- Je ans t ail of a BB and co r respo nd to 



a ~ 2. As clearly illustrated by iPanagia &: Fellil (119751 ) and IWright fe Barlowl (119751 ) (see 



also iPanagial 1 1 99 ll ) . ionized stellar winds display flatter SEDs than BBs. This is because at 
low optical depths the gas emission in the wind (mostly free-free emission) is rather flat, 
and in the optically thick regime the increase of the size of the emitting region as a function 
of wavelength (e.g. R oc A ' 7 for a constant velocity wind) partially compensates for the 
intrinsic thermal A~ 2 behaviour. Wind SED indexes asymptotically reach a = 0.6 in the 
case in which both the expansion velocity is constant (e.g. at suitably large radii where the 
terminal velocity has been attained) and the optical depth is high (e.g. long wavelengths). 

(iii) A model comprising a 30,000 K blackbody and an ionized wind. This approxi- 
mates the case of all early-type stars (essentially OB and WR types) whose winds add a 
flatter spectral component that is noticed as an excess which becomes m ore conspicuous at 



longer wavelengths. Calculations were made following IPanagial (119911 ) prescriptions for a 
spherically symmetric wind whose expansion velocity increases with radius like a power law, 
v = v a (r /ro) 1 ', and for a wide range of mass-loss rates. Although these power laws represent 
just a zero-order approximation to the actual velocity structure of the wind (most noticeably, 
lacking a limiting/terminal velocity behaviour as suggested by ultraviolet line studies), they 
are able to reproduce the power-law behaviour of the wind SEDs in the regime of high optical 
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depths, i.e. where the wind contribution becomes dominant. For the sake of simplicity, we 
show here only the curves corresponding to 7 = 0.5 and 7 = 0. Since the gas opacity is 
an increasing function of the wavelength, for each given mass-loss rate, the wind emission is 
higher at long wavelengths. As a result the corresponding trajectories on the various TCDs 
are curves that start at the colors appropriate for a pure stellar SED (a power law a ~ 2), 
end at the asymptotic slope for very optically thick winds (e.g., 0.95 for 7 = 0.5 and 0.6 for 
7 = 0), and display a characteristic upward concave curvature because of a more prompt 
increase in the redder colors (i.e. a positive second derivative). This is the only model not 
labeled in Figures [7] and due to space constraints. 

(iv) A model comprising a 30,000 K blackbody plus emission from an optically thin 
ionized nebula. Since the angular resolution of Spitzer is never better than l."7, which at 
the LMC distance corresponds to more than 0.4 pc, we have to consider the case of an 
unresolved HII region surr ounding an early-t y pe sta r. For these models we adopted the gas 



emissivity as computed by iNatta fc Panagial (119761 ). which includes emission lines, bound 



free and free-free continuum for both H and He, and is almost constant over the entire 
wavelength range 1 — 20 fim. As expected, the curves start from the bare atmosphere point 
and asymptotically tend to a point corresponding to a spectral slope of a ~ 0, displaying an 
upward concave curvature. 

(v) A model with a 30,000 K blackbody plus cool dust emission (140 K and 200 K BB) 
with the amount of dust increasing in logarithmic steps of 0.5. This model may illustrate the 
case of an early-type star surrounded by a cool dusty envelope, perhaps similar to a sgB[e] 



star (IKastner et al.ll2006l ). or stars with cold dust in their vicinity. In these calculations the 
dust opacity is assumed to be inversely proportional to the wavelength (re^ust oc A" 1 ). At 
these dust temperatures the peak of the dust emission occurs at about 18 and 26 /im, and, 
therefore, the corresponding dust SED is rapidly (almost exponentially) increasing toward 
longer wavelengths, thus displaying a very negative slope (a « 0). 

(vi) A model with a 3,500 K blackbody plus 250 K dust, which was selected to represent 
the case of red super giants (according to t he new RSG effective temperature scale for a 



M3 I in the LMC, see iLevesque et al.ll2007l ). Also in these calculations the dust opacity is 
assumed to be inversely proportional to the wavelength (ndust oc A -1 ), and the amount of 
dust is increased in logarithmic steps of 0.5. As expected, the starting point corresponds 
to the bare stellar atmosphere and the ending point is where the dust emission dominates 
(the peak of the dust emission SED occurs at about 15 /im), and the tracks display a typical 
upward concave curvature. 



From these simple models and considerations, it appears that ionized gas emission pro- 
duces a SED that declines with wavelength less steeply than a blackbody, thus possibly 
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causing an excess both in the infrared and in the radio domains but without reversing the 
overall declining trend of the SED. It follows that an observed flux increase toward longer 
wavelengths necessarily requires the presence of an additional component that has a suitably 
low temperature (e.g. less than 500 K for an upturn around 5 fim) and has a suitably high 
emission to affect, and eventually dominate the observed SED. These properties correspond 
to the characteristics of a dusty layer, that is able to absorb mostly optical and ultraviolet 
photons and reradiate that energy in the infrared. 

The locations of WR, LBV, sgB[e] and RSG are discussed in §6], d E] and [9J respectively. 



5. Infrared Excesses in OB stars 



The primary cause of mid-infrared excess in OB stars is free-free emission from winds 
or disks. Since this emission is a monotonically increasing function of the optical depth that 
depends on the square of the electron density, one may find a different mid-infrared emission 
from stars with identical ma ss-loss rates but diff erent wind velocity fields. Therefore, winds 
that are slow, cl umped fe.g. jBlomme et al.ll2003h or compressed toward the equatorial plane 
to form a disk (jOwocki et al.l 1 19941) will have s t ronge r mi d-infrared emi s sion. Given the 



sensitivity to these parameters, ICastor fc Simonl (119831 ) and lAbbott et al.l (119841 ) concluded 
that mid-infrared emission alone cannot be used to determine the mass-loss rates of OB 



stars. More recently, iPuls et al.l (120061 ) argued that this very sensitivity makes mid-infrared 
emission extremely important in disentangling the nature of clumped winds. We therefore 
proceed to examine the mid-infrared excesses of the 354 O and 586 early-B stars in our 
photometric catalog. 

In Figures 191 and fTOl we plot Jirsf vs. Jirsf — [3.6], Jirsf — [5.8] and Jirsf — [8.0] colors 
for O and early-B stars, respectively, indicating their luminosity classes, binarity and emission 
line classification properties by different symbols. We compare t he observed colors with colors 
of plane-parallel non-LTE TLUSTY stellar atmosphere models (ILanz fc Hubenyll2003l . 120071 ) 
of appropriate metallicity and effective temperatures. We have not dereddened the stars, 
because this requires BV photometry of higher accuracy than that provide d by the MCPS 



catal og; however, the low extinction toward the LMC {E(B — V) = 0.14 mag; iNikolaev et al. 



20041 1 makes such corrections small, as illustrated by the reddening vectors. For reference, 
TLUSTY models reddened by E(B — V) — 0.2 mag are also shown. We clearly detect 
infrared excesses despite not having dereddened the stars. At longer wavelengths, the excess 
is larger because the flux due to free-free emission for optically thin winds remains essentially 
constant with wavelength. Fewer stars are detected at longer wavelengths primarily because 
of the decreasing sensitivity of Spitzer and the overall decline of the SED. 
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We find that for J < 13 mag, more luminous OB stars exhibit larger infrared excesses 
in all colors due to winds. This can easily be understood considering that the infrared 
excess is an increasing function of the infrared optical depth, which in turn is proportional 

In early-type stars the mass-loss rate is 



(e.g. 



Panagialll991f ). 



to the ratio M 2 /(v 2 exp Rl) 
found to be proportional to a power k (higher than unity) of the stellar luminosity, i.e. 
M cx L k with k ~ 1.2 — 1.6 (ILamers fc Leithererl Il993l ; IScuderi et al.lll998l ). The average 
expansion velocity is a function of the effective temperature, rangin g from an initial value 
at th e photosphere v ~ v soun d cx T^ 2 to a terminal value cx T e 2 ff ( Panagia fc Macchetto 
19821 ). so that approximately v exp ~ (t>ot>oo) ' 5 oc T c g 25 . Since the stellar radius can be 



expressed as R cx L 1 / 2 /T e 2 ff , it follows that the infrared optical depth in early- type star winds 
is proportional to L( 2fc ~ 1,5 ) x T e ^ 5 , i.e. a ~ 0.9 — 1.7 power of the luminosity. This result can 
straightforwardly account for the observed increase of the infrared excess as a function of 
luminosity for O and B type stars in the LMC. We would like to mention that our findings 
confirm and pro yide a simple explanation for the luminosity effect between hypergiants and 
dwarfs noted by iLeitherer fc Wolfi (119841 ) in their sample of 82 Galactic OB(A) stars. The 
brightest stars with the largest excesses shown in Figure [9] are late O supergiants, with large 
M and low v^. The spread in excesses at any given J— band magnitude (~ 0.2 mag at 
J — [3.6], ~ 0.4 mag at J — [8.0]) cannot be due to reddening, which would displace all OB 
stars in a similar fashion. It is most likely related to the range of mechanisms that produce 
free-free emission, the properties of the winds (M, Vqo, clumping factors) and stellar rotation 
rates. 

Stars with J > 13 mag that exhibit large excesses in Figure [9] are classified as emission 
line Oe/On/One, i.e. with evidence of fast rotation or circumstellar disks. In Figure fTUl a 
larger fraction of stars with J > 13 mag are found to have large exce sses, reflecting the highe r 
fraction of emission line stars among B-stars compared to O-stars (INegueruela et al.l 120041 ) . 
The "B extr" stars ( Schild 1966 ; Garmany fc Humphreys 1985 ; Conti et al. Il986l) are found 
to occupy similar parts of the CMP as the Be(F e II) stars (notation by lMassey et al.lll995l ) 
and Be-Fe stars (notation by lEvans et al.l 120061 ) (and also Be/X-ray binaries), suggesting 
these are objects of one and the same nature, with the presence of Fe II lines perhaps due 
to higher densities in the circumstellar disks. We note that some Be stars do not exhibit 
excesses, while several stars not classified as Be stars have excesses. Variability (disk dissipa- 
tion), disk orientation, or insufficient spectral coverage leading to inaccurate classifications 
can explain these outliers. Therefore, infrared CMDs of optically blue selected stars can be 
used to identify new Be stars. An alternate explanation for some o f the outlier s is the new 
class of double periodic variables (or DPVs; iMennickent et al.ll2003l . 120051 . 120081 ). DPVs are 
evolved B and A-type giants in an Algol mass-transfer configuration with circumprimary and 
circumbinary disks. These disks are thought to be responsible for the observed long period 
variability and infrared excess detected in the near-infrared. 
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5.1. OB star SEDs 

Ultimately, the position of a star in a two-color or color-magnitude diagram is the result 
of its SED and luminosity, and it is the SED that can help us interpret its position in terms 
of its physical properties. Using the 0.3 — 24 /im photometry from Tabled we created SEDs 
for the stars by c onverting the i r mag nitudes to fl uxefl We used ef fective wavelengths and 



calibration s from iBessell et all (119981 ) for UBVI, iRieke et all (120081 ) for the 2MASS JHK» 
and IRSF, iReach et al.l J2005h for the Spitzer IRAC bands, and the MIPS Data Handbook 



for the MIPS 24 /im band. We compared the observed SEDs with models to assess the 
infrared excesses in the OB stars. TLUSTY models for LMC metallicity, a micro-turbulent 
velocity of 10 km s -1 , and the lowest available surface gravity for each effective temperature 
T eff were selected to correspond to supergiant stars. The exact value of log g has little effect 
on mid-infrared fluxes for models with T e R > 25kK and, while larger, the effect remains 
< 0.1 mag in (J — A) (for A > 2.2 /im) for the cooler models. We did not attempt to fit the 
SEDs, but simply overplotted models to provide a benchmark for determining the deviation 
of the observed SEDs from a "bare" photosphere. Deviations of OB stars from a TLUSTY 
model are due to reddening and free-free emission from their stellar winds. 

Figure [11] shows representative SEDs for 10 OB stars, normalized by the flux at the 
J— band. The MCPS, IRSF and SAGE measurements (filled circles), and the 2MASS and 
OGLE photometry (open circles) are connected by a solid line. TLUSTY models, similarly 
normalized, are overplotted as solid lines, while the dotted lines correspond to the models 
reddened by E(B — V) — 0.25 and 0.50 mag (see the Appendix for details on the treatment 
of interstellar extinction). In Figure [121 we plot the same SEDs, divided by the unreddened 
TLUSTY models. These unitless plots make deviations from models obvious by avoiding the 
enormous flux ranges encountered when plotting the SEDs of hot stars in the mid-infrared. 
This ease of inspection comes with a price: when plotted against logarithmic wavelength, 
equal bins on the abscissa do not correspond to equal energy intervals. However, in these 
plots, the effects of reddening (clearly seen in thet optical), infrared excess (detected in most 
cases) and even variability between the 2MASS-IRSF and OGLE III-MCPS photometry are 
clearly seen. 

The O stars shown in the left panels of Figures [11] and [12] include the lightly reddened, 
luminous and hot star Sk —70° 91, (02 III(f*)+OB), which only has a weak infrared ex- 
cess, in spite of having very strong ultraviolet P-Cygni wind lines (v^ = 3150 km s -1 ; e.g. 



6 Note, both magnitudes and fluxes exist in the SAGE database, but only magnitudes are included in our 
photometric catalog. Table 4 can be used to convert to fluxes. 

7 http:/ /ssc. spitzer. caltech.edu/mips/dh 
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Massa et al.ll2003l ). In such stars, the wind density remains low, creating only weak infrared 



emission. In contrast, the 04-5 V((f))pec star Sk —70° 60 and the 06 Iaf+ star Sk —65° 22 
(with = 1350 km s _1 ) have strong infrared excesses. The following two stars: Sk —66° 110 
and BI 192, are lightly reddened 09 I and 09 III stars, respectively, with distinctive evi- 
dence for infrared emission from their winds. These examples demonstrate that mass-loss 
rates cannot be derived solely from infrared excesses; an estimate of the spectral type and 
luminosity class is also required. We show examples of B giant and supergiant stars in the 
right panels of Figures [TT] and H2J The BO la star Sk —68° 52 is not consistent with the SED 
of a reddened 30kK model, but shows evidence for infrared emission from a wind. Contrast 
the BO III star Sk —67° 210, which appears to be nearly unreddened, to the moderately 
reddened BO III star ST2-46, exhibiting infrared excess. Sk —68° 8 (B5 Ia+), compared to 
a 15kK model, similarly shows an infrared excess. The final SED is of Sk —66° 12 (B8 la), 
also compared to a 15kK model (the coolest TLUSTY model available). While the star is 
probably cooler than 15kK, the reduction of the effective temperature by a few kK cannot ac- 
count for the strong infrared excess, which is consistent with free-free emission from a wind. 
Careful modeling of each star will be necessary to accurately derive their wind properties. 



6. Wolf-Rayet stars 



Of the 125 WR stars in the catalog of iBreysacher et al.l (119991 . not including BAT99-106 
and BAT99-108 through BAT99-112 in the crowded R136 cluster), 99 yield matches in the 
SAGE databasfl However, BAT99-45 a nd BAT99-83 are LB Vs and BAT99-6 is an O star 
binary (given its 03f*+0 spectral type; iNiemela et al.ll200ll ). Only 3 of the remaining 96 
WR stars were detected at 24 /mi: BAT99-22 (WN9h), BAT99-55 (WNllh), BAT99-133 
(WNllh, a spectroscopic twin to the LBV BAT99-83 in minimum). Our catalog includes 
various WR subtypes found in the LMC: early and late WN stars (some with hydrogen 
present); early WC stars of WC4 subtype; and one W03 star. While most WR stars are 
thought to be post-main se quence stars w ith progenitor masses greater than 25 — 75 M Q 
(depending on the subtype; ICrowtherl 120071 ). there is increasing evidence that the late-type 
WN stars with hydrogen in their spectra are hydrogen- burning stars with in itial masses 
above 40 — 60 (the lower initial mass be ing u ncertain: ISmith &: Conti 2008 ). Studies by 



Bartzakos et al.l (120011 ) , iFoellmi et al.l (120031 ) , and ISchnurr et al 



(120081 ) found the fraction of 



spectroscopic binaries among the WR stars in the LMC to be lower than expected. Given the 
diverse properties of objects encompassed under the 'WR" classification, one might expect 
the position of WR stars on a CMD or TCD to depend on their spectral types and possibly 
be influenced by binarity. 



We label the B0 I + WN3 star Sk -69° 194 |Massev et al.ll2000l) as an early-B star in our plots 
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We examine the position of the WR stars on the J vs. J — [3.6] CMD (Figure [TBI) as 
a function of their spectral type, subdividing them into: WN2-5, WN6-7, WN8, WN9-11, 
WC4 and W03 and labeling the known binaries. On average, the WN6-7 and WN9-11 
stars are more luminous and therefore brighter at J, when compared to their early WN2-5 
counterparts. In particular, the late WN9-11 stars containing hydrogen are the brightest at 
J. The se are the most luminous evolved WN stars, th ought to be quiescent states of massive 
LBVs (IBohannan fc Walborn!ll989l ; ISmith et al.lll994l ). Overall, the WR stars span 6 mag in 
J, and 2 mag in J — [3.6] color, with the brighter stars being "bluer", i.e. with steeper SEDs 
(see below). The brightest star is BAT99-22 (R84 or Sk -69° 79) with a WN9h spectr al type, 
which is known to be blended with a red supergiant (jHeydari-Malayeri et al.l 119971 ). while 
the reddest star is BAT99-38 (HD 34602 or Sk -67° 1 04), a triple system (WC4(+0?)+08I:; 
Moffat et al.lll990l ) with an associated ring nebula (IDopita et al.lll994[ ). which is probably 
responsible for the infrared excess. Dust is known to form around some colliding wind binary 
WC s tars (see references inlCrowtherll2007l ). and was recently detected in the vicinity of WN 
stars (IBarniske et al.l 120081 ) . The 24 jum detections around the WNllh stars BAT99-55 and 
BAT99-133, if confirmed, provide additional examples of WN stars with associated dust. 

In Figure [T4"[ we examine the colors of the WR stars in a J — K s vs. K s — [8.0] 
plot, follo wing lHadfield et al.l (|2007l) . The location of LMC WR stars does not agree with 
the region lHadfield et al.l (120071 ) defined by dashed lines to represent colors of Galactic WR 
stars. There are two reasons for this: 1) the Galactic sam ple includes many dusty WC9 stars 
that are intrinsically reddened (IGehrz fc Hackwelllll974l ). and 2) line-of-sight reddening has 
scattered the Galactic sample into the region these authors have defined. Metallicity is not 
expected to affect the colors of WR stars, in particular WC stars, whose winds consist of 
processed and therefore enriched material. The WR stars define a rather narrow linear trend, 
independent of spectral type, defined by: < J — K s < 0.7 mag and 0.4 < K s — [8.0] < 2.0 
mag, corresponding to power law spectra with indexes a = 1.5 to 0.5, i.e. optically thick 
winds with modest velocity gradients (e.g. |Panagialll99ll ). Single WR stars are expected to 
have flatter SEDs than binaries, which typically have main-sequence companions with steeper 
SEDs. The outliers in this plot are once again the heavily reddened WC4 triple system 
BAT99-38 and the blended star BAT99-22 (R84), whose infrared spectrum is dominated by 
the red supergiant. R84 is an outlier in all the CMDs, located among the RSG. On the 
[3.6] - [4.5] vs. [4.5] - [8.0] TCD (Figure HD the WR stars follow a similar linear trend. We 
note that the brightest WR stars in J are the bluest in the TCDs. 

In Figure [Jj)J we show SEDs for 7 WR stars with representativ e subtypes. The strong 
C IV and He II emission lines present in the W03 and WC4 stars (see lBreysacher et al.lll999l ) 
are responsible for the "kinks" in the SEDs of these stars. For 3 of these stars (BAT99-123, 
BAT99-61, BAT99-17) we show model fits derived from mainly optical and near-infrared 
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spectr oscopy using the spherical, non-LTE, line-blanketed CM FGEN code of iHillier &: Miller 



rll998h . For BAT99-123 (W03), we use the model presented in lCrowther et al.l J200oh which 
is derived from line profile fits to ultraviolet and optical spectroscopy. This model has 
parameters of T eff = 150 kK, log L/L Q = 5.3 , ««, = 4, 100 kms" 1 and M = 1 x 1O- 5 M 
yr -1 for a volume fi lling factor of 0.1. For BAT99-61 (WC4), we use the model presented 
bv ICrowther et al.l hoO$ ) which has T eff = 52 kK, log L/L Q = 5.68, = 3,200 kms" 1 
and M = 4 x 10" 5 M yr" 1 for a volume fillin g factor of 0.1. For BAT99-17 (WN4o), we 
use the model summarized in ICrowtherl (120061 ) which has T e g = 52 kK, log L/L Q = 5.40 , 
Voo = 1, 750 kms -1 and M = 1.4 x 10 -5 M Q yr _1 for a volume filling factor of 0.1. Overall, 
the model fits presented in Figure [15] show that the observed SEDs in the mid-infrared are 
flatter than the model continua, suggesting that either the model mass loss rates are too low 
or the winds are more highly clumped than the models predic t. We caution, ho wever, that 
this comparison is only for 3 stars. Generic CMFGEN models (ISmith et al.ll2002l ) are shown 
for the other 4 WR stars: a WN 85kK model for BAT99-1 (WN3b) and a WN 45kK model 
for the remaining WNL stars, all of which show excess above that predicted by the models. 
It is clearly desirable to investigate these apparent differences further by performing detailed 
line profile fits in the mid-infrared. 

Finally, in Table [5] we present photometry for 10 WR stars with counterparts (within 
20") at 70 /im and for 1 WR star with a counterpart at 160 fim. The star name, spectral type, 
magnitude and associated error at 70 /im and 160 /im, followed by the flux and associated 
error in these bands, are given. We caution that the angular resolution of MIPS is 18" and 
40", respectively, at these bands. However, we note that the sgB[e] and LBVs detected also 
show evidence for dust at the IRAC+MIPS24 bands. 



Luminous Blue Variables 



There are 6 confirmed LBVs (see review by Humphreys &: Davidsonlll994l ) in the LMC: 
S Dor, BAT99-83 or R127, R 71, R 110, BAT99-45, a nd R 85, al though another 5-6 have 
been suggested as candidates in the literature (see e.g., IWeisll2003h . The LBVs are not only 
among the most luminous sources at 3.6 /im (see Figure [2] and [3]) , with [3.6]-[4.5] colors 
similar to AGB stars and intermediate between RSG and sgB[e] stars, but also at 8.0 /im 
(Figure HJ) and 24 /im (Figure ED- In the TCDs (Figures M and CO) the LBVs a re located 
between the OB and WR stars, with the exception of R71 (ILennon et al.lll993l ). which is 
an extreme outlier, with K s — [8.0] = 3.4 mag and [4.5] — [8.0] = 2.5 mag. Its recent 
brightening of 1.5 mag between 2006 and 2009 (see light curve from the All Sky Automated 
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S^rue^ Pojmanski 2002 ) cannot account for its infra red colors. Instead , the emission from 
polycyclic aromatic hydrocarbons (PAHs) detected by lVoors et al.l (119991 ) in R71 can explain 
its "red" colors, since there are strong PAH lines in the [3.6] and [8.0] bands. All 6 confirmed 
LBVs were detected in the IRAC bands, 3 were detected in the MIPS 24 /im band and 2 at 
70 fim: R71 (saturated at 24 fim) and BAT99-83 or R127 (see Table EJ. The R71 de t ection 
is consistent with the 60 fim value in the IRAS Point Source Catalog and lVoors et al.l (Il999l ) 
suggest that a combination of crystalline silicates and cool dust are responsible for it. The 
detection of R127 at 70 fim may imply that a similar dusty environment exists around this 
star. We note that BAT99-45 is distinctively "bluer" than the other two LBVs detected at 
8.0 /iin and 24 fim (see Figures HI [5j |SJ). 

The SEDs of all 6 confirmed LBVs are shown in Figure [TS] and, unlike other classes of 
stars, are highly non-uniform. Since LBVs are by definition variable and the photometry 
in different bands was obtained at different epochs, the differences in the 2MASS and IRSF 
photometry are most likely due to variability. We only have MCPS photometry for two 
LBVs, which seems either variable or of suspect quality. S D or, the prototype S Doradus 
variable, and R127 (in outburst during the SAGE observations; IWalborn et al.ll2008l ) are the 
brightest LBVs at 24 fim (Figure E]). Many LBVs are surrounded by small nebulae (0.2-2 
pc) originating from past giant erupt i ons, which form cool dust. For example, R127 was 
found to have a bipolar nebula (jWeisI 120031 ) . The dust associated with its nebula is likely 
responsible for the bright mid-infrared magnitudes of R127. R71 exhibits a strong excess in 
the IRAC bands, as expected from its position in the TCDs. 

Do LBVs represent a homogeneous class of objects? Their infrared SEDs point to diverse 
properties, which is not surprising given their spectral variability during eruptions. On the 
one hand, the low luminosity LBV R71, with an estimated progenitor mass of around 40 M Q 
has one of the strongest infrared excesses of any massive star in the LMC and is well known 
as a source of dust emission. S Dor and R127 are more typical massive LBVs, both of which 
exhibit an infrared excess, possibly linked to circumstellar material formed during their recent 
outbursts. The other three LBVs in the sample: R110, BAT99-45 (also Sk -69° 142a, S83) 
and R85 all show different but coherent infrared behaviour more consistent with a stellar 
wind contribution only. It is possible that dust around these objects is cooler; however these 
stars were not detected at the longer MIPS bands. The different shapes of the LBV SEDs 
are likely related to the time since the last outburst event and the amount of dust formed. 



3 http: / /www. astrouw.edu.pl/asas/ 
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8. B[e] supergiants and Classical Be stars 



B[e] supergiants (sgB[e]) are a di stinct class of B-ty pe stars exhibiting forbidden emis- 
sion lines, or the "B[e] phenomenon" (ILamers et al.lll998l ). They are characterized by strong 
Balmer emission lines, narrow permitted and forbidden low-excitation emission lines and 
a strong mid-infrared excess. A two-component wind model with a hot polar wind and a 
slow and cool equatorial disk wind has been proposed to explain t hem (Zickgraf et all 1986): 



however models of the disk have so far achieved limited success ( Porter 



2003 



Kraus et al. 



20071 : IZsargo et al.ll2008l ). On the Hertz sprung- Russell (HR) diagram, sgB|e] stars are located 
below the Humphreys-Davidson limit ([Humphreys fc Davidson! Il979l ). a few of them being 
coincident with the location of LBVs. However, the existence of an evolutionary connec- 
tion of sgB [e| stars to LBV stars (transitional phase between Of and WR stars) is unclear. 
Langer fc Hegerl (119981 ) outline three possible ways sgB[e] could form circumstellar disks: 
from single massive evolved stars close to critical rotation (lMeynet &: Maederl 120061 ). blue 
supergiants which have le ft the red supergiant branch, and from massive binary mergers 
(IPodsiadlowski et al.ll2006l ). 



The sgB[e] stars are the most conspicuous group of stars in all infrared CMDs and 
TCDs: they are amon g the brightest and most reddened stars in the LMC (also noted by 
Buchanan et al.l I2006T ) . In the LMC, 12 stars have been c lassified as sgB [e] stars (including 
LH 85-10, although it is not among the 11 stars listed in IZickgrafl 120061 ) . The 11 that were 
included in our catalog (S 22, S 134, R 126, R 66, R 82, S 12, LH 85-10, S 35, S 59, S 137, 
S 93) were all matched in the SAGE database. Their SEDs are shown in Figure [161 The 
B0.5[e] stars are sorted by decreasing flux at [3.6], the later-type B[e] stars by decreasing 
effective temperature. A 25kK TLUSTY model representing the underlying photosphere is 
overplot ted. We find t hat th e SED of LH 85-10, a suspected LBV assig ned a B|e] spectra l 
type by lMassey et al.l (120001 ) . resembles the SEDs of "classical Be" stars (INegueruelal 120041 ) . 
Figures [21 El El and [7| demonstrate that LH 85-10 is much fainter and less reddened than 
sgB[e] and LBV stars, and is located near Be stars, and therefore may have been misclassified 
as a B[e] star. The SEDs of the remaining 10 sgB[e] stars are all very similar, with slowly 
decreasing flux in the optical, an inflexion point in the near-infrared and a "bump" starting 
at 2 fim and peaking near 5 /im. This peak corresponds to hot dust at ~600 K. The slight 
change in the slopes of the SEDs between 8 and 24 /mi from star to star suggests different 
contributions from cool dust (150 K). R66, S35 and R126 were also detected at 70 /zm 
(Table EJ. We find all 10 sgB[e] stars to be extincted by E(B - V) > 0.5. 



Kastner et al.l (120061 ) presented Spitzer IRS spectroscopy for 2 of these sgB[e] (R66, 
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R12g 10 l). finding evidence for massive circumstellar disks, which can simultaneously explain 
the (relatively) small amount of reddening in the optical and the observed infrared excess. 
However, the mass-loss rate required to sustain such a massive disk is very high and unlikely 
to persist during the He burning phase of the star. These authors argued against a model 
of episodic mass ejection based on the similarity of the SEDs of the two sgB[e] they studied. 
The similarity of all 10 sgB[e] SEDs makes the scenario of episodic mass ejection even less 
likely or, alternatively, suggests a shorter duration of the sgB[e] phase or even an origin in a 
binary companion. The small fraction of sgB[e] stars among B stars mig ht present a clue to 



their origin. We note that the rare, dust obscured, luminous progenitors (IPrietq et al 



2008 



Prietol 120081 ) of the new ly discovered class of transients, inc luding SN2008S (jSmith et al. 
2009|) and NGC 300 OT (Bond et alil2009l ; feerger et al]l2009T i. could be related to the lower 



luminosity B[e] supergiants, given that both are likely post-main- sequence stars (possibly 
in a post-red supergiant phase) with dust and similar luminosities. iThompson et al.l (120081 ) 
presented evidence for a short duration (< 10 4 yr) of the dust-enshrouded phase prior to 
eruption. 

The infrared emission in the sgB[e] stars is much stronger than in classical Be stars 
(shown in Figure [T7|) . It is typically several hundred times stronger than the photosphere at 
8 /im, and up to 1000 times stronger at 24 /im, explaining the extreme positions of sgB[e] 
stars in the TCDs. In the classical Be stars, free-free emission from disks can cause the 
8 fim flux to be 10 times that of the underlying photosphere, as is the case for all but 
one (Sk —71° 13) of th e Be stars shown. The e mission in classical Be stars can be highly 
variable (see review by iPorter fc Riviniud 120031 ) . which may account for the disagreement 
between the IRSF and 2MASS photometry for NGC 2004-033 (B1.5e), and in some cases, 
can even dissipate altogether, as for Sk —71° 13 (Be (Fe-Be)). The SEDs of the remaining 
Be, "B extr" , Be (Fe-Be) stars are identical, thereby providing additional evidence (see §5]) 
that these stars are the same objects: the most luminous of the Be stars. We also present 
the SED of 1 of the 4 Be/X-ray binaries in our matched catalog, finding it to have a peculiar 
shape, similar to the other Be stars, but with excess emission in the [3.6] and [4.5] bands, 
possibly due to variability. 



10 Note, the sgB[e] star R126 (Sk -69° 216) is located in NGC 2050, close to the LBV R127 (Sk -69° 
220), which is in NGC 2055. These two clusters could be a "fossil" two-stage starburst. 
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9. Red Supergiants & Other Stars 



The red K and M supergiants are among the brightest stars in the LMC at all infrared 
wavelengths. However at 24 /im they span 5.5 magnitudes (see Figure EJ). This is due to 
the different types and amounts of dust contributing to the SEP at t h ese w avelengths due 



to different mass-loss rates (see iVerhoelst et al.ll2009l ). IJosselin et al.l (120001 ) note that the 
K — [12] color index is a mass-loss indicator; therefore, we similarly interpret the range in 
[24] magnitudes and the 2 magnitude spread in K s — [8.0] color in Figure [6] as a reflection 
of the range of mass-loss rates. In Figure [HJ we notice a temperature trend: red supergiants 
clustered at the bottom left (i.e. with "blue" colors) are predominantly K or early M super- 
giants, while the bulk of M0-M4 supergiants are found in a band, roughly following the RSG 
model described in £JD The temperature trend correlates with the optical depth of the dust 
emission and the amount of dust, which in turn is a function of the mass- loss. The offset 
from the blackbody model can be explained by dust; exceptions to this trend might be due 
to misclassifications, blends, or variability. The same temperature trend is seen in Figure 
with the coolest RSG being the most luminous at 24 /im. Our catalog includes LMC 170452 , 
which is a rare example of a RSG that has changed spectral types. iLevesque et al.l (120071 ) 
found it to vary from M4.5-M5 I in 2001 and 2004 to Ml. 5 I in 2005 December. The SAGE 
observations were also taken in 2005 (July & October), possibly during the transition. We 
suggest that supergiants not appearing to follow the temperature trend in Figure M could 
signify spectral variability (e.g. the "blue" M2 star [SP77] 55-20 and the very red K0 star 
165543). The extremely reddened (at long wavelengths) outlier in Figures HI [5], and [8] is the 
M2-3 supergiant 175549 with [8] - [24] = 4 mag. The M0 I 138475 has K s - [8.0] = 3 mag, 
making it an outlier in Figure [61 In Figure [3X it is located with the AGB stars, implying it 
is an M-type AGB star. 



Figure [18] shows representative SEDs of yellow and red supergiants. The SED peak 
moves redward with decreasing temperature, as expected. Note that starting at Sk —69° 30 
(G5 I), a depression in the continuum due to the CO band appears at 4.5 /im, and persists 
to later types. This depression is extremely weak or absent in the Ml la, 138552. Some of 
the later supergiants have a distinct break in their slopes at roughly 4.5 /im, also coinciding 
with the CO band. The SED of L H 31-1002 (F2 I ), sug gests a hotter effective temperature, 
inconsistent with its spectral type. iMassey et al.l (120001 ) give B — V = 0.4 mag for this star, 
in disagreement with the value from MCPS (B — V — —0.2 mag), perhaps suggesting an 
error in its coordinates that lead to a match with a nearby hot star. The K7 I star 139027 
also appears to have an exce ss at the shortest wavelengths; therefore it is a candidate ( 
Aurigae binary (jWrightl 1 19 701 ). with a hot companion. LMC 170452, has a peculiar SED 
(with V — K = 7.2 mag, instead of ~ 5), possibly reflecting its spectral type change. 
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Finally, in Figure [19] we present some peculiar SEDs, illustrating examples of nebular 
contamination (Nll-081), variability (RXJ0544. 1-7100), blending (BAT99-22) or misidenti- 
fication (Sk —67° 29, Sk —70° 33) in our photometric catalog. Nebular contamination and 
misidentifications are responsible for the OB outliers in the CMDs and TCDs. Only 3 out 
of 354 O stars (Nl 1-028, Nl 1-046, lmc2-703- all main-sequence stars) and 15 out of 586 
early-B stars have 24 /zm detections, illustrating that contamination is less than 3% at this 
wavelength. 



10. Summary 

This paper presents the first major catalogs of accurate spectral types and multi- 
wavelength photometry of massive stars in the LMC. The spectroscopic catalog contains 1750 
massive stars, with accurate positions and spectral types compiled from the literature, and in- 
cludes the largest modern spectroscopic surveys of hot massive stars in the LMC. The photo- 
metric catalog comprises uniform 0.3—24 fim photometry in the UBVIJH i^ s +IRAC+MIPS24 
bands for a subset of 1268 stars that were matched in the SAGE database. Our photometric 
catalog increases the infrared photometry for OB stars presented in a single study by an order 
of magnitude, allowing for a detailed study of their wind parameters. The low foreground 
reddening toward the LMC and the identical distance of the stars in the sample remove 
degeneracies inherent in Galactic studies and enable the investigation of infrared excesses, 
while minimizing systematic errors due to reddening. 

We examine the infrared excesses of the half-solar metallicity OB stars in the LMC, 
finding a correlation of excess with wavelength, as expected, and luminosity. The dispersion 
in the relation reflects the range of mechanisms, and therefore stellar wind parameters, which 
produce free-free emission in OB stars. We discuss the positions of OB stars, WR, LBV, 
sgB[e], classical Be stars, RSG, AFG supergiants and Be/X-ray binaries on CMDs and TCDs, 
which serve as a roadmap for interpreting existing infrared photometry of massive stars in 
nearby galaxies. We conclude that sgB[e], RSG and LBVs are among the most luminous 
stars in the LMC at all infrared wavelengths. Representative SEDs are shown for all types of 
stars, and in certain cases are compared to model atmospheres to illustrate infrared excesses 
due to free-free emission from winds or disks. We confirm the presence of dust around 
10 sgB[e] stars from the shape of their spectral energy distributions, which are presented 
for the first time at these mid-infrared wavelengths, and find the SED shapes to be very 
similar. The large luminosities and the presence of dust characterizing both B[e] supergiants 
(log L/L Q > 4) and the rare, dusty progenitors of the new class of optical transients (e.g. 
SN 2008S and NGC 300 OT), suggest a common origin for these objects. The variety of 
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SED shapes observed among the LBVs is likely related to the time since the last outburst 
event and the amount of dust formed. Finally, we find the distribution of infrared colors for 
WR stars to differ from that of Galactic WR stars, due to the abundance of Galactic dusty 
WC9 types and the effects of foreground reddening in the Milky Way 

This work demonstrates the wealth of information contained in the SAGE survey and 
aims to motivate more detailed studies of the massive stars in the LMC. Examples of studies 
that can be conducted using our catalogs follow: (a) Candidate RSG, sgB[e], LBV in the 
LMC and other galaxies can be identified by selecting them from infrared CMDs, ruling out 
AGB stars from optical and near- infrared photometry (if it exists), and confirming these 
with spectroscopy. Fo r example, in IC 1613, a comparison of the IRAC CMD in Figure 4 of 
Jackson et al.l (l2007bl ) to our Figure [H suggests that the brigh t est "b lue" stars are candidate 
RSG or LBV stars. With the optical catalog of iGarcia et al.l (120091 ). their visual colors can 
be determined and 0.3-8 /zm SEDs can be constructed, (b) Specific stars or categories of 
stars can be selected from the catalog for further study. For example, mass-loss rates can be 
deter mined for OB stars t hat have spectra from the Far Ultraviolet Spectroscopic Explorer 



(e.g. iFullerton et al.ll2000l ). since terminal velocities can be measured from the ultraviolet 
lines, (c) The OB stars can be used to determine the LMC reddening law, given the known 
spectral types and SED shapes of these stars, (d) Variability in 2MASS-IRSF photometry 
can be used to pick out new, variable, massive stars. The most luminous stars in our catalog 
(V ~ 10 mag) are saturated in the LMC microlensing surveys; however, the two epochs of 
near-infrared photomet ry can be used, for examp le, to find ca ndidate eclipsing bin aries not 
included in the OGLE JWvrzvkowski et al.l 120031 ) or MACHO iFaccioli et alil2007h lists, (e) 
Our catalog can be cross correlated with X-ray catalogs of the LMC to study non-thermal 
emission processes in massive stars. 

A comparison of the infrared properties of massive stars in the Large vs. Small Magel 



lanic Cloud (SMC) will be pursued next. Stellar winds scale with metallicity (IMokiem et al. 



20071 ); therefore, a comparison of infrared excesses of OB stars in the LMC vs. the SMC will 
help quantify the effect of metallicity on OB star winds, Be star properties, sgB[e] stars etc. 
In the future, the James Webb Space Telescope will obtain infrared photometry of massive 
stars in galaxies out to 10 Mpc and measuring the colors of massive stars as a function of 
metallicity, will be useful for interpreting these data. 
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Appendix 

We adopted a reddening curve composed of a continuous contribution and a feature 
which represents sil icate absorption at 9.7 nm . The continuous portion is based on the 
formulation used by lFitzpatrick &: Massal (120091 ) to fit near- infrared extinction cu rves (whic h 
include 2MASS JHK S ) and is a generalization of the analytic formula given by iPeil ( 119921 ). 
For A > Ao, it has the form 



k(X-V) c = 



E(X - V) 0.349 + 2M7R(V) 



E(B - V) 



(A/Ac 



- my), 



where A = 0.507 /im. The ratio of total to selective extinction, R(V) = A(V)/E(B — V), 
and a are free parameters. Near-infrared observations show that these parameters can 
vary from one sight line to the next. In this work, we use a curve defined by a = 2.05 and 
R( V) = 3.11. These parameter s provide a good representation of the mean curve determined 
by iFitzpatrick fc Massal (120091 ) and also to extinction data taken from the literature for 
diffuse interstellar medium sight lines in the Milky Way . The latter data s ets extend to 
24 /im and include the fo l lowin g: ICastor fc Simonl (119831). iKoornneel Jl983j ) jAbbott et c 
Jl984h - iLeitherer & Word Jl984h . iRieke et all Jl989h . IWeeneri Jl994h . and iNishivama eTal. 



(120091 ). Although we cannot be certain this same form applies to the LMC, it is consistent 



with the SEDs of lightly reddened, normal hot stars that are thought to have weak winds (see 
§2]). In any case, the extinction for most of our objects is small enough that inaccuracies in 
the interpretatio n of reddening should be m inimal. For completeness, we included a Drude 



profile (see, e.g., IFitzpatrick fc Massal 120091 ) to describe the 9.7 /im silicate feature, even 



though it has little effect on the Spitzer bands. The adopted profile is given by 



/(A) 



\x' 



x 



2\2 



+ (x'yf 



(2) 



where x = l/X /im 1 , and the parameters were assigned the values xq = 1/9.7 /im 1 , 7 = 
0.03 and 03 = 3 x 10~ 4 . Finally, the complete extinction curve is given by k(X — V) — k(X — 
V) c + f(X). When we plot the SEDs, we actually use the form k(X-J) = k(X-V) -k( J-V). 
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Fig. 1. — Spatial distribution of massive stars with IRAC counterparts, which are found to 
trace the spiral features in the 8 /im SAGE mosaic image. Different symbols denote different 
spectral types. The spatial distribution of the OB stars is uniform, with the exception of 
the Nil (at ~ 74°, -66.5°), NGC 2004 (at - 82°, -67.3°) and 30 Doradus (at ~ 85°, -69°) 
regions. 
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[3.6]-[4.5] 



Fig. 2. — [3.6] vs. [3.6] — [4.5] color magnitude diagram for massive stars with IRAC coun- 
terparts in the SAGE database. The c onversion to absol ute magnitudes is based on a true 
LMC distance modulus of 18.41 mag (iMacri et al.l 120061 ). Different symbols denote differ- 
ent spectral types. The locations of all the SAGE detections are shown in grey as a Hess 
diagram. The sgB[e], RSG and LBVs are among the most luminous stars at 3.6 /im. 
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J-[3.6] 

Fig. 3. — Same as FigureEJ but for the [3.6] — J vs. [3.6] color magnitude diagram. The longer 
baseline separates the populations more clearly. The reddening vector for E(B — V) = 0.2 
mag is shown. 
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[8.0]-[24] 



Fig. 4. — Same as Figure HI but for the [8.0] vs. [8.0] — [24] color magnitude diagram. The 
sgB[e], RSG and LBVs are also among the most luminous stars at 8 jum. 
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[8.0]-[24] 



Fig. 5.- 
brightest 



Same as Figure [21 but for the [24] vs. [8.0] — [24] color magnitude diagram. The 
sgB[e], RSG and LBVs are among the most luminous stars at 24 /im. 
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Fig. 6. — J — K s vs. K s — [8.0] diagram for massive stars in our catalog. The locations 
of all the SAGE detections are shown in grey as a Hess diagram. The solid lines represent 
models (described in §1]): (i) a BB at various temperatures, as labelled, (ii) a power law 
model F„ oc v a , for —1.5 < a < 2, (hi) an OB star plus an ionized wind (not labelled), 
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[g-*]-[9-e] 

Fig. 7. — Same as Figure El but for the [3.6] — [4.5] vs. [4.5] — [8.0] diagram. The majority 
of hot massive stars lie between the blackbody and OB star +wind model, illustrating that 
a BB is a good approximation in the infrared. 
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Fig. 8. — Same as Figure El but for the [3.6] — [8.0] vs. [8.0] — [24] diagram. The RSG with 
"bluer" colors have earlier spectral types, indicating a temperature sequence (see $9]). 
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Fig. 9. — Infrared excesses {Jirsf vs. Jirsf — [3-6], Jirsf ~ [5-8] and Jirsf — [8-0]) for 354 
stars. Supergiants are shown in yellow, giants in green, main-sequence stars in blue, stars 
with uncertain classifications ( "other" ) in red, binaries with a large circle and Oe stars with 
an x. The solid lines correspond to 30kK and 50kK TLUSTY models with logg = 4.0. A 
reddening vector for E(B — V) — 0.2 mag is shown, as well as reddened TLUSTY models by 
this same amount (dotted lines). The more luminous stars exhibit larger infrared excesses, 
which increase with A. The 0.2 — 0.4 mag spread in excesses at any Jirsf— band magnitude 
reflects the range in mass-loss rates, terminal velocities, clumping properties and, perhaps, 
rotation rates of O stars. 
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Fig. 10. — Same as Figure [9j but for 586 early-B stars. The solid lines correspond to 20kK, 
hgg = 3.0 and 30kK, log# = 4.0 TLUSTY models. A reddening vector for E(B - V) = 0.2 
mag is shown, as well as reddened TLUSTY models by this same amount (dotted lines). 
The larger number of early-B stars makes the trends identified among the O stars clearer. 
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Fig. 11. — Representative SEDs of stars (Left) and B stars (Right), normalized by their 
J— band fluxes (dashed line) and offset for display purposes. Normalized TLUSTY model 
atmospheres (30kK, log g = 3.0 model for the O and BO stars; 15kK, logg = 1.75 model for 
the late B stars) are overplotted for comparison. The MCPS, IRSF and SAGE measurements 
are shown as filled circles; the 2MASS and OGLE measurements as open circles. The dotted 
curves correspond to TLUSTY models reddened by E(B — V) — 0.25 and 0.50 mag. Infrared 
excesses are detected in most stars (see §5. 1 [) . 
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Fig. 12. — Ratios of the SEDs of O stars (Left) and B stars (Right) to the unreddened 
TLUSTY models shown in Figure [TT], clearly showing the deviations of the observations from 
the models and cases of variability, from a comparison of VIJHK S magnitudes obtained from 
different sources. Dotted curves correspond to TLUSTY models reddened by E(B — V) = 
0.25 and 0.50 mag. Infrared excesses are present in most cases. 
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Fig. 13. — J vs. J — [3.6] color magnitude diagram for Wolf-Rayet stars with infrared 
detections, labelled according to spectral type. Black circles denote known binaries. The 
reddening vector for E(B — V) = 0.2 mag is shown. The WN2-5 stars are on average fainter, 
while the WN6-7 and WN9-11 are brighter at J— band. 
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Fig. 14. — J — K s vs. K s — [8.0] diagram for Wolf-Rayet stars with infrared detections, 
showing a tight linear c orrelation indep e ndent of spectral type. The dashed lines correspond 
to the colors found by lHadfield et al.l (120071 ) for Galactic WR stars; the solid line shows 
a power law model oc u a , for < a < 2. The difference in colors between the LMC 
WR stars and the Galactic WR stars is a combination of extinction, as illustrated by the 
reddening vectors, and the large fraction of dusty WC9 stars among Galactic WR stars. 
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Fig. 15. — Same as Figure [TTl but for Wolf-Rayet stars (Left) and Luminous Blue Variables 
(Right). For the WR stars, we show model fits for BAT99-123, BAT99-61 and BAT99-17 
and generic CMFGEN models for the rest (see $6]). All show excess above that predicted 
by the models. The first 3 LBVs show evidence for dust; BAT99-83 and R71 (saturated at 
[24]) were detected in the MIPS70 band. The different shapes of the LBV SEDs are likely 
related to the time since the last outburst event and the amount of dust formed. 
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Fig. 16.— Same as Figured! but for B[e] supergiants, with a 20kK, logg = 2.25 TLUSTY 
model overplotted. The sgB[e] have similar SEDs revealing dust, except for LH 85-10, which 
has a SED similar to that of a Be star. R126, R66, and S35 were also detected in the MIPS70 
band (see §SJ). 
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Fig. 17. — Same as Figure [TT1 but for classical Be stars, with a 25kK, log g = 2.75 TLUSTY 
model overplotted. The similarity of their SEDs implies that the various spectral types refer 
to the same type of object. The Be/X-ray binary RXJ0535. 0-6700 exhibits excess emission 
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Fig. 18. — Same as Figure [TTJ but for A, F and G supergiants (Left) and K and M su- 
pergiants (Right). No reference models are shown. The SED of LH 31-1002 (F2 I) implies 
a hot effective temperature, inconsistent with its spectral type. The K7 I star 139027 has 
an excess at the shortest wavelengths, possibly indicating a hot companion. The M4.5-5 
supergiant 170452 is a rare case of a RSG that changed spectral type (see 
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Peculiar Stars 
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Fig. 19. — Same as Figure (TTJ but for peculiar SEDs, arising from nebular contamination 
(Nll-081), variability (RXJ0544. 1-7100), blending (BAT99-22) or misidentification at one 
(Sk -67° 29; at 24 /im), or all (Sk -70° 33) bands. 
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Table 1. CATALOG OF SPECTRAL TYPES FOR 1750 LMC MASSIVE STARS 



Star 


RA (J2000) 


Dec (J2000) 


Reference 6 


Classification 


Name" 


(deg) 


(deg) 




& Comments 


BAT99-1 


71.3843765 


-70.2530289 


F03 


WN3b Sk -70 1 


Sk -69 1 


71.3865433 


-69.4670029 


C86 


Bl III 


Sk -67 1 


71.4227524 


-67.5705261 


C86 


B extr 


Sk -67 2 


71.7685394 


-67.1147461 


F88 


Bl Ia+ 


Sk -70 la 


71.8649979 


-70.5626984 


C86 


09 II 


Sk -69 8 


72.3293762 


-69.4318314 


F88 


B5 la 


BAT99-2 


72.4009000 


-69.3484000 


F03 


WN2b(h) 


Sk -69 9 


72.4642105 


-69.2011948 


J01 


06.5 III 



a Star des ignations: Erevsacher etaD (BAT99: Il999h ISanduleaH ( Sk: Il970h . 
Westerlundl (W; ll96ll) . lBrunet et al.1 (Bl: Il975h , iLuckel (LH;[l972j), iHenize! fS: 

1956) 



Re ference key : (B0 1) iBartzakos et all (1200 lh . (B09) iBonanosI ( 12009) 



(C86) IConti et all dl986l). (E06)lEvans et al 



1991), (F02) Fitzpatrick et al. (120021) . (F03a) 



Fitzpatrick et all (120031) . (F03) 



20061) . (F88) i iFitzpatrickl d 19881 



Foellmi et al 




(G98) 



(200 31), (G94) 



Guinan et al 



(199 



Garmanv et all (119941). (O05) 



Gonzalez et al 



(H92) Hcvdari-Malav eri fe Melnick 



. (H03) Hcvdari-Malaveri et al 
Humphreys fc Davidson! dl994l) . (J01) 



( 20031). (H79)lHumphrevsldl979l). (H94) 



d2005l). (M95 ) iMassev et al.1 dl995l) 



Jaxon et all (I2001D. CL05) 



(M00 ) 
( 20021) 



Massevl (120021) (M0 2b) IMassev et al 
Massev et all (I2004I). (M05) IMassev et al 



Niemela et al 



(OOlb) lOstrov et all (12 001) 



2003) 
19931) 



(P01) 
(R02) 



(Me05) iMevnadier et all (12005^ (N94) 



Massev et all (|2000D (M02) 



Liu et al 



(2003J), (M04) 



(M03) iMasse 
( 20051). (M091 IMassev! (|2009l ). 
Niemela fc Bassinol dl994) (N01 ) 



(2001). (095) lOev fc Massev 



'arker et al 



Schild fc Testor 



2001), (P92) 



1995), (O01) 



Olsenetal 



(OOl c) lOstrovl (120011). (O03)lQstrov fc Lapasset 



(200 21), (R05) 



'arker et al 



(2001) 



(1992), (P93) iParker 



liba set al 

dl992l). ( S 08) ISchnurr et al.l (I2008D. (T98) 



Raeuzova fc Popovl (120051). (S92) 



1998 ), (W77) 



Walborn et al 



2002 ), (W03) 



Testor fc Niemela 



Walbornl dl977l). (W97) IWalborn fc Blades! dl997l). (W02) 



Walborn et al 



(2004), (W08) IWilliams et al.l ( 2008h . (Z06) lzickgral (|2006h 



(2003) , (W0 4) IWalborn et al 



Note. — Table Q] is available in its entirety in the electronic version of the 
Journal. A portion is shown here for guidance regarding its form and content. 
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Table 2. Statistics for the 1268 Matched Stars 



Catalogs Matched # Stars 



IRACC 5 

IRACC+IRSF 88 

IRACC+IRSF+MCPS 601 

IRACC+IRSF+MCPS+MIPS24 122 
IRACC+IRSF+MCPS+MIPS24+OGLE 9 

IRACC+IRSF+MCPS+OGLE 364 

IRACC+IRSF+MIPS24 21 

IRACC+IRSF+OGLE 18 

IRACC+MCPS 33 

IRACC+MCPS+MIPS24 1 

IRACA+MIPS24 1 

IRACA+MCPS+MIPS24 1 

IRACCEP1+MCPS+MIPS24 4 



Table 3. 0.3-24 jjm Photometry of 1268 Massive Stars in the LMC 



StarName Q IRAC Designation RA(J2000) Dec(J2000) U a v B a B ... Ref. b Classification & Comments 



BAT99-1 


J044532.25 


-701510.8 


71.3843765 


-70.2530289 


Sk -69 1 


J044532.76 


-692801.2 


71.3865433 


-69.4670029 


Sk -67 1 


J044541.44 


-673413.9 


71.4227524 


-67.5705261 


Sk -67 2 


J044704.44 


-670653.1 


71.7685394 


-67.1147461 


Sk -70 la 


J044727.58 


-703345.6 


71.8649979 


-70.5626984 


Sk -69 8 


J044919.03 


-692554.6 


72.3293762 


-69.4318314 


BAT99-2 


J044936.24 


-692054.8 


72.4009 


-69.3484 



15.001 


0.037 


15.647 


0.029 .. 


. F03 


WN3b Sk -70 1 


12.343 


0.01 


13.364 


0.01 


. C86 


Bl III 


12.829 


0.038 


13.72 


0.025 .. 


. C86 


B extr 


10.351 


0.104 


11.219 


0.095 .. 


. F88 


Bl Ia+ 


12.693 


0.01 


13.654 


0.01 


. C86 


09 II 


10.838 


0.06 


11.714 


0.078 .. 


. F88 


B5 la 


15.633 


0.038 


16.495 


0.03 


. F03 


WN2b(h) 



a Star de signa tions: iBrevsacher et all (BAT99; 1 19991) . ISanduleakl (Sk; Il970h . IWesterlundl (W; Il96ll) . iBrunet et aD (Bl; Il975l ). iLuckel (LH; Il9721 ). 
Henizel CS: ll956l) 

b Re ferenc e key: ( 1301) iBartzakos et all d200ll) . CB0 9) iBonanosI d2009l ) (C8 6) IConti et all dl986l) (E06) lEvans et all d2006l). CF881 IFitzpatrick! 



Gonzalez et al 



Massev et al 



( 20051) . ( G98) 



Guinan et al 



(1998) 



(1988, 19911). (F 02) IFitzpatrick et all (120 02). (F03a) iFitzpatrick et all (120031). CF03) iFoellmi et all (120031) CG94) iGarmanv et al 



Humphreys (1979), (H94) Humphreys & Davidson (1994), (J01) 



2000), ( M02) IMassevI (20021). (M02b) 



Massev et al 



19951) COOl) lOlsen et al l (l200lh. /OOlb) lOstrov et all (120011). fOOl c 



Schnurr et al 



Walborn et all (|2003h . (W04) 



(2005 
d2001 



(H92) IHevdari-Malaveri fc Melnickl (119921). fH03) iHevdari-Malaveri et al 



Walborn et al 



>ey 

(|2005|) . ('M09) lMassevl (I2009D. (M e05) iMevnadier et all (l2005h . CN94) iNiemela fc BassinT(ll994D. (NOT) 



Jaxon et all (120011). (L 5) iLiu et all (120051). CM95) iMassev et al 



(2002) . (M03) IMassevI (120031) . fM0 4) IMassev et all (l2004h . (M05) 



Nicmela et al 



CP92) barker et all il9m . /P93) IParkerl jl993T). fR 02) iRibas et all J2OO2I) fR05) iRaguzova fc : Popovl (hooi) 



(Il994h . fG6jg) 
(2003), (H79) 
dl995l ). (MOO) 



)lQstrovl d20 0lh. rO03)lQstrov fc LaoassetJ d2003h. fPOl) 



(120081). ('T98) lTestor fc Niemelad 19981) . CW77 ) IWalbornl dl977h CW97) I Walborn fc Blades! (fl997h . fW02) IWalborn et all d2002I) . (W03) 



Massev et al 



(|200ll ). r095)l0ev fc Massev 



d2()()lli. fW08) IWilliams et all t>l)()A). CZ06) IZickgraj (|2006l) 



Parker et al. (|200lh . 
CS92) ISchild fc Testorl dl992l) . (S08) 



Note, 
content. 



Table [3] is available in its entirety in the electronic version of the Journal. A portion is shown here for guidance regarding its form and 
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Table 4. Filter & Detection Characteristics 



Filter 


Acs 


Zero mag 


Resolution 


# stars 




H 


flux (Jy) 


(") 


detected 


U 


0.36 


1790 


1.5/2.6 


1121 


D 


0.44 


4063 


1.5/2.6 


1136 


V 


0.555 


3636 


1.5/2.6 


1136 


I 


0.79 


2416 


1.5/2.6 


1026 


VoGLE 


0.555 


3636 


1.2 


391 


loGLE 


0.79 


2416 


1.2 


391 


J 


1.235 


1594 


2.5 


1203 


H 


1.662 


1024 


2.5 


1218 


K s 


2.159 


666.7 


2.5 


1184 


JlRSF 


1.235 


1594 


1.3 


1122 


HlRSF 


1.662 


1024 


1.3 


1089 


Ks IRS f 


2.159 


666.7 


1.3 


1077 


[3.6] 


3.55 


280.9 


1.7 


1260 


[4.5] 


4.493 


179.7 


1.7 


1234 


[5.8] 


5.731 


115.0 


1.9 


950 


[8.0] 


7.872 


64.13 


2 


577 


[24] 


23.68 


7.15 


6 


159 



Table 5. MIPS70+160 Photometry 



Star Name 


Spectral Type 


mag7o 




FIUX70 


O F/ 1/T 7n 




Flux 1 Rn 


(7 fin 






(mas;) 

V 0/ 


(mag) 


(mjy) 


(mjy) 


(mae;) (mag) 

\ 0/ V 0/ 


(mjy) 


(mjy) 


BAT99-8 


WC4 


1.379 


0.02796 


218.4 


5.623 








BAT99-22 


WN9h (+RSG) 


2.327 


0.04657 


91.24 


3.911 








BAT99-32 


WN6(h) 


1.4 


0.03326 


214.4 


6.564 








BAT99-37 


WN3o 


1.503 


0.02847 


194.8 


5.107 








BAT99-38 


WC4+[08L] 


-0.8254 


0.01853 


1664 


28.4 


-2.412 0.1344 


1476 


181.8 


BAT99-53 


WC4 (+OB) 


1.386 


0.03978 


217.1 


7.952 








BAT99-55 


WNllh 


0.6481 


0.03829 


428.3 


15.1 








BAT99-85 


WC4 (+OB) 


-0.4681 


0.03374 


1197 


37.19 








BAT99-123 


W03 


1.782 


0.04259 


150.7 


5.909 








BAT99-133 


WNllh 


0.3187 


0.01856 


580.1 


9.914 








R66 


B0.5[e] 


0.4201 


0.01875 


528.4 


9.126 








S35 


Bl[e] lab 


2.104 


0.05634 


112 


5.806 








R126 


B0.5[e] 


0.993 


0.03794 


311.7 


10.89 








R71 


LBV 


-0.9245 


0.01292 


1823 


21.7 








R127 


LBV 


-0.1097 


0.03093 


860.7 


24.51 









